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Abstract. We explore the formation and evolution of debris ejected around quark stars in the Quark Nova scenario, and the 
application to Soft Gamma-ray Repeaters (SGRs) and Anomolous X-ray Pulsars (AXPs). If an isolated neutron star explodes 
as a Quark Nova, an iron-rich shell of degenerate matter forms from its crust. This model can account for many of the observed 
features of SGRs and AXPs such as: (i) the two types of bursts (giant and regular); (ii) the spin-up and spin-down episodes 
during and following the bursts with associated increases in P\ (iii) the energetics of the boxing day burst, SGR 1806+20; (iv) 
the presence of an iron line as observed in SGR1900+14; (v) the correlation between the far-infrared and the X-ray fluxes 
during the bursting episode and the quiescent phase; (vi) the hard X-ray component observed in SGRs during the giant bursts, 
and (vii) the discrepancy between the ages of SGRs/AXPs and their supernova remnants. We also find a natural evolutionary 
relationship between SGRs and AXPs in our model which predicts that the youngest SGRs/AXPs are the most likely to exhibit 
strong bursting. Many features of X-ray Dim Isolated Neutron stars (XDINs) are also accounted for in our model such as, (i) 
the two-component blackbody spectra; (ii) the absorption lines around 300 eV; and (iii) the excess optical emission. 
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1. Introduction 

Soft y-ray Repeaters (SGRs) are sources of recurrent, short 
(f ~ 0.1 s), intense (L ~ 10'*'* ergs) bursts of y-ray emission 
with a soft energy spectrum. The normal patterns of SGRs are 
intense activity periods which can last weeks or months, sep- 
arated by quiescent phases lasting years or decades. The three 
most intense SGR bursts ever recorded were the 5 March 1979 
giant flare of SGR 0526-66 (Mazets et al. 1979), the similar 28 
August 1998 giant flare of SGR 1900-H14 and the 27 December 
2004 burst (SGR 1806-20). AXPs are similar in nature but 
with a somewhat weaker intensities and no recurrent bursting. 
Several SGRs/AXPs have been found to be X-ray pulsars with 
unusually high spin-down rates of P/P ~ 10"'" s"', usually at- 
tributed to magnetic braking caused by a super-strong magnetic 
field. 

Occasionally SGRs enter into active episodes producing 
many short X-ray bursts; extremely rarely (about once per 50 
years per source), SGRs emit giant flares, events with total 
energies at least 1000 times higher than their typical bursts. 
Current theory explains this energy release as the result of 
a catastrophic reconfiguration of a magnetar's magnetic field. 
Magnetars are neutron stars whose X-ray emission are pow- 
ered by ultrastrong magnetic fields, B ~ 10'^ G. Although the 
magnetar model has had successes, in this paper we present an 
alternative model which addresses outstanding questions. 



We explore these issues within the quark-nova (QN) sce- 
nario (Ouyed et al. 2002; Ouyed et al. 2004; Keranen, Ouyed, 
& Jaikumar 2005). In our previous studies we have suggested 
that CFL (color-flavor locked) quark stars could be responsi- 
ble for the activity of SGRs and AXPs (Ouyed et al. 2005; 
Niebergal et al. 2006), and in this paper we extend the QN 
model by studying in more details the evolution of the QN 
ejecta (as first discussed in Keranen & Ouyed 2003). 

Following the QN explosion we show that a high metal- 
icity shell can form from the neutron star crusQ. If this mat- 
ter gains sufficient angular momentum from the central quark 
star it can form a torus-like structure via the propeller mech- 
anism, which is discussed in the second paper in this series 
(Ouyed et al. 2006; hereafter refered to as Paper II). However 
if it does not, as we discuss in this paper, there is instead a thin 
corotating shell suspended by the quark star's magnetic field 
pressure. Because the quark star is in a superconducting state 
it's magnetic field decay is coupled to its period in a manner 
prescribed by Niebergal et al. (2006), such that as the field de- 
cays the shell will drift closer to the star. This movement will 
shift sections of the shell above the line of neutrality, loosely 
defined as the polar angle (measured from the equator) above 
which the magnetic field vector is sufficiently parallel to the 
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' The words "suspended crust" has been used in a different context 
already with respect to quark stars; i.e. the electro-magnetically sus- 
pended crust at hundreds of Fermi's above the quark star surface (e.g. 
Alcock et al. 1986) instead of a few quark star radii as in our case. 
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gravity vector to allow material to break off from the shell, and 
fall into the star. Upon collision with the star the pieces of the 
shell are instantly converted to CFL quark matter (described by 
Lugones & Horvath (2002)), where excess energy is released 
as high-energy radiation (Ouyed et al. 2004). We propose this 
as the mechanism responsible for SGR and AXP bursting ac- 
tivity, and also show that SGRs and AXPs differ primarily by 
age. Another class of objects are also explored, namely X-ray 
Dim Isolated Neutron stars (XDINs; e.g. Haberl 2004), and we 
find that these may have evolved from SGRs/ AXPs. 

This paper is presented as follows. In § |2] we review the 
concepts of a quark-nova and the expansion/evolution of the 
ejected material into a shell. § [3] contains discussions on the 
geometry of the resulting shell and its self-similar behaviour 
in time, leading to pieces of the shell breaking off. The sub- 
sequent high-energy bursting from these pieces falling into the 
quark star is then described in § [H along with methods for es- 
timating the age of SGRs/AXPs. In §|5]the quiescent phase of 
SGRs/ AXPs is discussed in the context of our model, where a 
simple relation between luminosity and spin-down rate is de- 
rived. The changes in period and period derivative observed 
during bursts are then discussed within the framework of our 
model in § |6l § [7] contains case studies for specific SGRs 
(1806-20 and 1900H-14), along with explanations for the pres- 
ence of an iron line during bursts, correlated X-ray to infrared 
flux ratio, and the hard spectrum seen in giant flares. We then 
summarize in § 17.51 some outstanding issues in the current un- 
derstanding of SGRs/AXPs, and show how our model can pro- 
vide explanations. Finally in §[8] we discuss XDINs and show 
how many features are readily explained using our model. We 
then conclude in §|9l 

2. Quark Nova 

In the quark-nova (QN) scenario, the core of a neutron star 
(NS) shrinks to a stable, compact, quark object before the con- 
version of the entire star to (u,d,s) matter. By contracting, and 
physically separating from the overlaying material (hadronic 
envelope which is mostly made of crust material), the core 
drives the collapse (free-fall) of the overlaying matter leading 
to both gravitational and phase transition energy release as high 
as 10^^ ergs in the form of neutrinos. The result is a quark star 
in the Color-Flavor Locked (CFL) superconducting phase, sur- 
rounded by an ejected shell. Although it has been shown that 
pure CFL matter is rigorously electrically neutral (Raj ag opal & 
Wilczek 2001), other work (Usov 2004 and references therein) 
indicates that a thin crust is allowed around a quark star due 
to surface depletion of strange quarks. In our model we have 
assumed no depletion of strange quarks, which implies a bare 
quark star. 

In the CFL phase there are a total of nine combinations 
between the gluons and the photon with eight of these combi- 
nations subject to the Meissner effect. The only one that does 
not suffer the Meissner effect involves a combination of elec- 
tromagnetism and a U(l) subgroup of the color interactions 
(e.g. Ferrer et al. 2005). In other words the magnetic field 
will penetrate some phases of the color field but not others. 
Unfortunately, the relevant calculations are done in effective 



models of QCD which cannot accurately handle the details of 
the mixing and the corresponding back reaction of the quarks. 
Assumptions had to be made (e.g. Meissner effect) in order to 
proceed with our astrophysical model which implies that we 
are ignoring the one component that penetrates the supercon- 
ductor. In the CFL quark star, given our basic assumptions, the 
magnetic field is contained in the vortex array so the internal 
field is uniform and equal to the surface field. 

The QN ejecta consists mainly of the neutron star's metal- 
rich outer layers heated by neutrino bursts. It was shown that 
up to IO^^Mq can be released during the QN (Keranen, Ouyed, 
& Jaikumar 2005). Thus our model proposes a high density, 
metal-rich ejecta surrounding a superconducting quark star. 

Within our model there are four possible scenarios as to the 
outcome of the ejecta. First, if the ejecta is very light then it will 
become gravitationally unbound from the quark star, where 
the r-process begins creating heavier elements (Jaikumar et al. 
2006). Second, if the ejecta is too heavy it will fall back into 
the quark star releasing tremendous amounts of energy. This 
we propose could lead to explosive events reminiscent of short 
gamma-ray bursts, which we are currently exploring. Third, is 
the case where the ejecta mass is such that it can be suspended 
above the surface by the quark star's magnetic pressure. Fourth, 
the ejecta was formed with enough angular momentum to move 
into a Keplerian orbit. The third case is discussed in this paper, 
whereas the fourth is discussed in the second paper of this se- 
ries. 

2.1. Expansion of the ejecta 

We consider the dynamics of the shell (of mass msh) ejected 
during the QN. If the shell's outward momentum is not too 
large, the interplay between gravity and the QS's magnetic 
field is dominant. Thus, as the highly conducting shell expands 
spherically outwards, an electromotive (/ x B) force is induced 
to oppose its motion, causing it to continuously decelerat^ 
The natural oscillatory motion created by the magnetic pres- 
sure gradient and gravity is resultingly damped out, resulting 
in the ejecta finding an equilibrium radius where the forces due 
to the magnetic pressure gradient and gravity are in balance. 
This equilibrium radius is then given by, 

^ lOkmA/sinflp ■ (1) 

-y/msh, -6^^05,1.4 

Here we have used - Bs(^.s/^m)^ for a dipole magnetic 
field with a strength at the surface of, Bs.is^ in units of 10'^ G. 
The star's mass and radius are in units of 1.4Mq and 10 km 
respectively, while the mass of the shell, msh.-6, is in units of 
IO^^Mq. The line of neutrality, 0^, is measured from the equa- 
tor and defines the angle above which the gravity vector is no 
longer mostly perpendicular to the magnetic field vector If sec- 
tions of the shell are above this line, they are free to break off 

^ The damping term is by far the dominant term in the equations of 
motion of the shell due to the high conductivity of the shell, estimated 
by 0" = iicshe^A/ (in^Cs). Here «e,th 's the number density of thermal 
electrons, Cj is the sound speed, the mean free path is given by A = 
1/ ("e,thO"T)> and (Tj is the Thompson scattering cross section. 
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and fall into the star's poles along the field lines. The condition 
for the shell to exist is that it is off the surface of the quark star, 
> ^Qs, which imphes ^^js^qsjo ^i°^B > (OT.sh,-6A^Qs,i.4)- 
For this equilibrium situation to be achieved the mass and 
initial velocity of the shell are confined to a somewhat nar- 
row range of values. Although the mass of the shell is self- 
consistent with previous studies of neutron star crust masses 
(discussed below), the initial velocity is not similarly con- 
strained. As such, we assume it can take on a range of values, 
with the small initial velocity case corresponding to the work 
presented in this paper, and the high velocity case is left for 
the second paper in this series. We feel that this dichotomy is 
physical, and explains some observed AXP features well. 

Recent work studying natural mechanism of magnetic field 
amplification in quark matter (just before or during the onset 
of superconductivity; Iwazaki 2005), shows that 10^^ G mag- 
netic fields are readily achievable. This amplification can occur 
in quark matter due to the response of quarks to the sponta- 
neous magnetization of the gluons. In contrast, for the magne- 
tar model, the sole proposed mechanism for generating 10'^ G 
fields requires millisecond proto-neutron stars. This seems to 
be challenged by recent observations (Allen & Hovarth 2006). 

Thus, equation ([T]i implies that for magnetic fields in the 
range of 10'"^-10'^ G, the maximum mass of the shell cannot 
exceed lO^^-lO"*" Mq. A higher shell mass translates into < 
Rqs, meaning that the magnetic field is not sufficiently strong 
to stop the shell from promptly falling back onto the star This 
situation is very possible but not considered in this paper, as it 
is left for future work. 

The shell's size is determined by the NS crust density 
profile (Datta et al. 1995), and is of the order of (0.012- 
0.025)] X Rqs corresponding to shell masses of the order of 
(lO-'*-lO-*)M0 for densities of (1.8 x 10^-5.5 x 10'°) g/cc. We 
note that if degenerate, the shell will always be in relativistic 
degeneracy since the densities are above pcr = 7.3 x lOVc g/cc 
(Lang 1974). Here, fi^ is the mean mass per electron and is 
taken to be //g ~ 2 since the shell's maximum density is below 
the neutron drip density, ~ 4 x 10" g cm"^ . 

The transition from Fermi-Dirac to Boltzmann statistics oc- 
curs at the degeneracy temperature Tfenni - 76.2 MeV pjj, for a 

1 /3 

non-relativistic gas and Tfevmi — 8.8 MeV pjj, for a relativistic 
gas (e.g. Lang, p 253). So, for the shell to be born degenerate 
at 10 MeV, it implies a minimum density of ~ lO'" g/cc which 
translates to an ejected mass of ~ 10"^ Mq. We note that even 
if the shell density at birth is less than ~ 10^° g/cc, (blackbody) 
cooling during the expansion of the ejecta to R^n results in a 
shell temperature of less than 1 MeV, thus, yielding a degener- 
ate shell. 

Once the shell reaches the magnetic radius it will have 
expanded to a thickness of, AR^, such that it satisfies hydro- 
static equilibrium. Using the relativistic degenerate equation- 
of-state (P - xp'^l'^ where x - 1.244 x lO'^yUg'*^"'), the width of 
the shell, positioned at ^m, can be calculated to be. 



A^ 



= 2.1 X 10 



m 



M, 



-3/4 



.shell,-6'"QS,1.4 



(2) 



This implies that the shell's thickness does not change much 
when it reaches Rm- The corresponding average density of the 
shell at the magnetic radius is then, 

,9 



,9/4 j^9/4 



p,„ = 7.58 X 10 g/cc m^j, _g.v.Q5 ^^x., ^5XXQs 10 
implying that the shell is relativistic degenerate 

2.2. Energy release 



(3) 



The internal energy of the expanding ejecta is U - 
N(3/5)eF 1 + 7f (^)^ where = ?Msheii//^e'WH, is the total 
number of electrons. The Fermi energy in this case is ep ~ 
8.8 MeWp\l^ (e.g. Pathria pg.200). As a resuh of the heating 
of the electrons, the optically thick shell will radiate approxi- 
mately as a blackbody. 



sh 



(4) 



The total energy it releases as it relocates to R^n is then (taking 
kT « ep) 



Af/p ~ 5.2 X 10^5 erg '«sh,-6(p;/io - P.'/jo) 



(5) 



where p, lo and pm,m are the densities of the shell initially and at 
magnetic equilibrium radius respectively, in units of lO'*' g/cc. 
This energy released during expansion is done so in roughly a 
millisecond. 

We have not investigated the fate of this energy release, 
whether it is absorbed by adiabatic expansion losses or radi- 
ated. What is crucial to our model is the fact that the shell 
remains degenerate at R^ which is guaranteed by hydrostatic 
equilibrium. 

2.3. Angular momentum transfer 

For a magnetic radius, R^, larger than the corotation radius 
the propeller will take effect (Schwartzman 1970; Illarionov 
& Sunyaev 1975), deflecting the QN shell into a torus on the 
equatorial plane. Using an angular momentum conservation ar- 
gument, we can estimate the location of such a torus by writing 



(6) 



where R^ is the equatorial location of the torus in a Keplerian or- 
bit around the quark star, and Qk = ^GMqs /R^ is the Kepler 
rotation. Applying equation OJ implies a torus radius of. 



p4 q2 
"m"QS 

Rt ^ 0.2 km 

GMqs 



sin^0B<i5^QS.lQ 
<-6^QS,1.4^?™s 



(7) 



where the birth period, Pi ms, is given in units of milliseconds. 

However, in order to actually form a torus we require 
enough angular momentum transfer to guarantee > ^m- This 
translates into an upper Umit on the initial period of, 



Aiim < 0.46 ms 



rsin fl ^3/4 d3/2 p9/2 
(Sm6yB) ^,j5^QS,i() 

sh,-6 QS,1.4 



(8) 



This presents two cases: we study the first in this first paper 
when the star's initial period is too large for torus formation. 
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This is where the shell will remain at and be in corotation 
with the magnetosphere. The second case where the shell can 
achieve a Keplerian orbit and form a torus is studied in paper II 
(Ouyed et al. 2006). 



3. Co-rotating Degenerate Shell 

3.1. Shell geometry 

The area of a shell at the magnetic equilibrium radius, is 
Ash = 47r/?^ sin 6b, where is the polar angle (measured from 
the equator) and defines the line of neutrality. In other words, 
for 9 > 6b the shell material is free to "slip" along the magnetic 
field lines onto the star, whereas for 6 < 6b the shell mate- 
rial is suspended by magnetic pressure. Thus, the geometry is 
such that there is a thin shell at the equator subtending an an- 
gle of 26b, and empty regions at the poles. The shell is shaped 
in the star's distorted magnetic dipole with an inward bulge at 
the equator, and the radius of the shell at the equator, under 
hydrostatic balance, can be shown to be roughly half that at 6b- 



3.2. Self-similarity: shell dynamics and evolution 

The radius of the shell, R^, is proportional to the magnetic 
field strength (i.e. Eq. [1]), so we use the magnetic field decay 
prescribed by Niebergal et al. (2006) to calculate the radius in 
time. This is because the quark star itself is in the superconduct- 
ing phase, implying that vortex expulsion couples field decay 
and period evolution. 

After a characteristic time, (as defined by Eq. 4 in 
Niebergal et al. 2006), the shell will respond to the decaying 
field by moving towards the star, resulting in a new equilib- 
rium radius. During this movement, sections of the shell will 
be shifted above the line of neutrality, 6b, where they will 
be broken off and fall into the star along the magnetic field 
lines. Upon colliding with the quark star these shell pieces will 
be converted immediately to CFL quark matter, releasing ex- 
cess energy in the form of radiation, which in our model is an 
SGR/AXP burst. We then have a new state, n + 1, where the 
shell has moved towards the star, lost mass, and is now sitting 
at radius ^m,n+i with a mass msh.n+i- 

The new characteristic decay time, period of rotation, pe- 
riod derivative, and magnetic field are respectively. 



Tn+l = K- 



p2 



1/3 



= ^(1 + 



5„+i = B„il + —) 



-1/6 



(9) 

(10) 

(11) 
(12) 



where the constant /<• = 8.8 x 10^** ■ s"', is for a quark star 
with a mass of I.AMq and radius of 10 km. Also, t, is the time 
elapsed since the previous burst (or birth of the quark star if it 



is the case). Taking the ratio of the new characteristic time to 
the old gives. 



I^ = l + - = 1 



(13) 



where a = t/T„ defines the typical decay time of the magnetic 
field for each epoch n. Because of the self-similar nature of 
the problem, a is roughly constant from one epoch to another, 
which implies that the time interval between two successive 
bursts. 



Af, 



aT„+i ^a(l+ a)"T„ 



(14) 



increases in time on average. 

So after a time t = aT„, the magnetic field decays and the 
shell originally sitting at R^ ^ moves in closer to the star Now 
by keeping the mass of the shell constant as it slowly drifts in 
from R^ to an inner radius Rin and applying equation ^ we 
can write. 



-1/6 



(15) 



In the above equation the notation used implies that the shell 
reaches Rin when the surface magnetic field has decayed to 
Bn+i- Following this, the shell will then suddenly lose mass 
at Ri„, due to the pieces of the shell closest to the poles being 
shifted above the line of neutrality, 6b, and breaking off. The 
resulting change of mass of the shell can be realized by con- 
sidering the area of shell material below the line of neutrality 
while keeping density in the shell roughly constant during the 
process. 



Wsh,n+i _ 47rsm6BRl _ 
rrisKn 47rsin0B^m,n ^m,« 



= (l+a)- 



1/3 



(16) 



The shell will now be less massive causing it to rapidly 
(i.e. keeping B constant, Bn+i = ^in) equilibriate by moving 
further away from the star to a new radius that defines the (n + l) 
state. It so happens, this new radius is the same as the original. 
This is seen from equation ([T]l, 



R 



m,n+ 1 



R„ 



B„ 



»^,sh,n 
'Wsh,n+1 



1 



(17) 



Therefore as the shell adjusts over time, it stays at roughly the 
original radius R^.n- In other words, it oscillates about the same 
position while losing mass. This mass loss can be determined 
appealing to equations (fT5]&[T6l). 



ATOsh,n = 'Wsh.n - '«sh,n+l = / (a) X OTsh.n 

where we define /(a) as, 
f(a)^[l-(l+ar''') . 



(18) 



(19) 



Thus, the amount of mass lost during a burst in this self-similar 
model (Eq. [T8] ) is determined by a. The corresponding energy 
release is. 



= 1.8x10'*^ erg 77o.i»J,sh,-6/(Q'), 



(20) 
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where 770.1 is the energy conversion efficiency factor taken to 
be ~ 0.1. We also note that the ratio between succesive burst 
energies is, 



n+l 



'«sh,n+I 



= (l+a) 



-1/3 



which implies that the bursts weaken in time. This is expected 
since the shell's mass is decreasing in time. 

In our model, AXPs are merely older versions of SGRs, 
thus one would expect SGR bursts to be more intense and fre- 
quent. As can be seen in table [Tl (where the objects are in or- 
der of decreasing estimated age), there is seemingly a decreas- 
ing trend in burst intensity and frequency as one moves from 
younger to older objects. It may be that both types of bursts 
have this trend, however, we do not feel there is enough events 
to draw any conclusions. This is a potential test for our model. 



4. The bursting phase 

The magnetic energy due to field decay is released continu- 
ously over a long timescale (thousands of years) and gives the 
steady x-ray luminosity of SGRs and AXPs (see |5). What is 
important and unique about the shell (i.e. the accretion energy) 
is that the energy is released in bursts, and the shell sitting at 
Rm offers a natural mechanism/torque for sudden changes in 
period derivative. 

After time t ~ ar, the magnetic field decays substantially 
and the entire shell moves in to Rin/Rm ~ (1 + a)"''*". A larger 
a means the shell moves closer to the star, causing larger sec- 
tions to be shifted above the line of neutrality, implying larger 
burstfl 

Another possible scenario is that giant bursts are due to 
global instabilities like Rayleigh-Taylor while regular bursts 
are due to chunks breaking-off the edge of the shell. We as- 
sume that the Rayleigh-Taylor instability does not act due to 
the long timescale for magnetic penetration of the conducting 
shell (see §4. 1 in paper 11). What is presented below is intended 
to present the overall scenario and follow-up work (i.e. numer- 
ical simulations) is necessary to investigate the effects global 
instabilities might have on this simplified picture. 



where /0.2 = f(a =1) and msh,n,-6 is the current mass of the 
shell in units of IO^^Mq. The waiting time between each bursts 
is on average. 



(21) Af = QT Si Q' X 2500 yr 



/ p2 



(23) 



where f 10 is the spin period in units of 10 s, and the magnetic 
field is in units of lO'^ G. One can see here that, because B 
decreases and P increases in time, the waiting times are in- 
creasingly less frequent as the object ages. Given the range in 
B derived for SGRs in our model (see Table 1), this implies a 
burst frequency of ^ l/(QflOOO yr) per object. 

4.2. Regular bursts: a «: 1 

If pieces of the shell are able to break off as the shell moves 
in by only a small amount, then a <K 1. This corresponds to 
small pieces of the shell falling into the star and, in our model, 
leads to regular-sized bursts. This process can be interpreted as 
small oscillations around R^ consistent with small fractional 
mass-loss episodes by the shell. Observationaly, regular bursts 
are known to follow a power-law distribution in energies which 
in our model would be due to a power-law distribution in a. 
Clearly, the detailed mechanism which determines a would be 
different for giant bursts and regular bursts, and would require 
numerical simulations to understand. 

For Of <K 1, equation ( fTSl l approximates to Awsh,n - 
{a/3)msh,n which yields burst energies of. 



~ a X (6 X lO"** erg s ') m.\m^,n,-b 



(24) 



Forff = 10"^, a burst energy of £1, ~ 6x10"*° erg is attained, and 
one should expect a waiting time of Af ~ 1 day (see Eq. l23T l. 
for a star with a 7 s rotation period and field strength of 10'^ G. 

Thus, in our model the difference between regular and gi- 
ant bursts is due to the size of piece that breaks off the shell. 
Although the size is of a stochastic nature, we expect younger 
objects (i.e. the SGRs) to have a more unstable shell, resulting 
in more large pieces breaking off, causing more frequent and 
intense bursts. For older objects (i.e. the AXPs) the shell be- 
comes more stable with age resulting in weaker, less frequent, 
bursts. 



4.1. Giant bursts: a ^ I 

If the time needed for the shell to move in sufficiently to have 
pieces broken off is roughly equal to the charactertic field decay 
time, then a = 1, and ~ 20% of the shell's mass is lost during 
the burst. According to equation |20l the corresponding energy 
release is, 

Eb,n = 3.6 X lO"**" erg 770.1 OTsh,n,-6/o.2 , (22) 

^ If Rin < Rqs then areas of the shell are in contact with the star. In 
this case, it is easy to imagine that the shell will experience a major 
disruption as inner sections are converted to CFL matter during con- 
tact probably destroying the entire shell; this could have applications 
to other explosive phenomena and is beyond the self-similar picture 
presented here. 



4.3. The P-P diagram 

From equations (fT0l&[T2b. the magnetic field is related to the 
period and period derivative by. 



B = Vi^. 



(25) 



Furthermore, combining this and the equation for (Eq. [TJ 
with equation ( l20b one can show that. 



sin 08 770.1/(0') ' 



(26) 



Illustrated in the upper panel of Figure [T] are cases with large 
and small a. The solid curve is for E\, - 10'^^ erg s^^, a - 2, 
and R^ = 20 km, and fits well with the boxing day event. Also 
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Fig. 1. The upper panel is a P-P diagram with SGRs/AXPs and XDINSs plotted. The solid curve is a plot of Eq. ( |26] | with a 
burst energy of - lO'** erg, a -2, and - 20 km. Also shown are curves for =100 km, a - 2 and energies Zib = 10'*'* 
erg s"' (dotted) and - lO'*' erg s"' (dashed). These line of constant energy move leftwards as an object ages. In the lower 
panel luminosity is plotted against period derivative according to Eq.[30l The upper and lower solid lines represent a magnetic to 
X-ray conversion efficiency, ?7x, of 1 and 0.01 respectively, and the dashed line is the luminosity averaged over all viewing angles 
(for //x = 1)- Since P cc 1/age, objects further to the left are generally older 



shown are two curves for R^^ = 100 km, a - 2 and energies 
£b = lO'*'* erg s-i (dotted) and £b = 10"*' erg s"' (dashed). The 
lines of constant burst energies in the upper panel of Figure [1] 
move to the left as time increases (as the object ages). This can 
be seen from equation {TH in § 13.21 which shows the burst en- 
ergy decreasing in time. We also note that as objects evolve fur- 



ther to the left, the bursts are weaker and less frequent, which 
is why AXPs in our model evolve from SGRs. 

4.4. Age estimate 

To determine the exact age of AXPs/SGRs in our model, one 
would need to know the entire bursting history of the object 
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which is impossible. For a real situation the burst sizes are vari- 
able implying a different a for each burst. So to approximate 
the age in the simplest case (when a is constant over time on 
average), the age elapsed since the QN is, 

n 

tQN -'^OTi . (27) 

i=I 

Using equation dTlT l the series can be summed to yield 

fQN=r„(l-(l+a)-") . (28) 

Provided that there has been a large number of bursts (n » 1), 

T„ = PJ3P„ (29) 

is a good estimate of the age, however we do expect a to vary 
from one burst to another. Although for SGRs, the energy in 
the large bursts {a ~ 1) dominates the small bursts (a ~ 10"^), 
resulting in an age very close to that of equation (|28] |. This age 
also happens to be the age estimate derived from equations (|9] 
- fTTT i when assuming f » r, which is reasonble for AXPs in 
our model. So our model predicts the age of ~ Q.5P/3P for 
SGRs, and ~ P/3P for AXPs and XDINSs. However, without 
knowing the burst history we cannot derive the initial period 
and magnetic field (see equations [TOl and fT2]) . 

One may be tempted to use the age of the associated parent 
supernova remnant to estime the age of the SGR/AXP, but in 
our model these ages are neccesarily different. This is due to 
the time required for a neutron star to reach quark deconfine- 
ment densities (Staff et al. 2006), as well as the time needed for 
strange quark nucleation (i.e. u,d to u,d,s) to occur (Bombaci 
et al. 2004). Together these delays can easily add up to the age 
difference between SGRs/ AXPs and their supernova remnants. 

5. The quiescent phase 

The X-ray luminosity during the quiescent phase of 
SGRs/ AXPs in our model is due to vortex expulsion from spin- 
down. The magnetic field contained within the vortices is also 
expelled, and the subsequent magnetic reconnection leads to 
the production of X-rays. From equation (22) in Ouyed et al. 
(2004) and equations (11011121 1. the luminosty emitted from this 
process is, 

Lx ^ 2.01 X 10^^ ergs-'7/xP^i, , (30) 

where rix is the efficiency parameter inherent in the conversion 
from magnetic energy to radiation. The lower panel in Figure[T] 
shows our model of luminosty versus period derivative as com- 
pared to the quiescent X-ray luminosities of AXPs, SGRs and 
XDINSs. 

The scatter in the AXPs/SGRs data can be explained within 
our model by considerations of shell geometry and viewing an- 
gle. More specifically, any X-rays being emitted from the star's 
surface towards the solid angle filled by the shell (47rsin6'B) 
will be affected due to absorption by the shell. X-rays are 
then able to escape from an open area at the poles with a 
solid angle on the order of (47r(l - sin 0b))- Because we sta- 
tistically observe viewing angles equally from all directions. 



this would imply an average (geometrical mean) lumimosity 
of Lx X (1 - sin^e) ~ 2.7 x IQ^'^rixP-ii ergs"', assuming 
= 60°. This is plotted by the dotted line in the lower panel 
of Figure [U Thus, the scatter in data is a combination of X-ray 
production efficiency and viewing angle effects. The two "out- 
liers" namely, AXP 1E2259H-586 and AXP 4U0142H-615 can 
be explained as those AXPs/SGRs born in the propeller regime 
so that they are surrounded by a degenerate torus instead of a 
shell, and are discussed in more detail in paper II. 

5.1. The two-component blackbody 

Our model provides an explanation to the two component 
blackbody spectrum that may be seen in SGRs/ AXPs (i.e. Israel 
2002). This is realized by considering the blackbody tem- 
perature as set by the X-ray continuum luminosity, Lx - 

4^^Qs7^BB't0b^ 

/ p \l/2 

rBB-0.41keV4/'* — ^ , (31) 

\"QS,10/ 

where we made use of equation (|30] | for Lx and adopted the 
fiducial value 7?Qs = 10 km. 

The effective shell temperature in our model is also 
set by the X-ray continuum luminosity, Lx sin % - 
sin Obo-T^^^^^ , and is 

/ p \i/2 

rsh,eff. - 0.41 keV 77]/' — . (32) 

\"m,l{)/ 

Thus there is two temperatures from two different emitting 
blackbodies, namely the quark star and the surrounding shell. 

Whether or not the blackbody component from the shell 
could actually be observed or not hasn't been considered. If it 
can not, then the magnetic field strengths around CFL quark 
stars is still at a magnitude that can produce the synchrotron 
break (i.e. blackbody plus power law) that is usually assumed 
for SGRs/AXPs. We only assert that if the blackbody emission 
from the shell could be observed, then its temperature is such 
that it correlates well with a two-component blackbody spec- 
trum model. 

6. Period changes and torques 

Angular momentum exchange between the shell and the star 
occurs during the shell's inward and outward movements, and 
would manifest itself as changes in rotational period and period 
derivative during bursts. We start with the inward drift where 
the moment of inertia of the star decreases, causing it to spin- 
up. 

6.1. Inward Shell Drift: Spin-up 

In our model the inward shell drift occurs while the vortex (i.e. 
magnetic flux) expulsion mechanism, as discussed in Niebergal 

■* Equation ( 131b is a direct measure of the star's radius once Tbb (or 
Lx) and P are measured. This could become crucial for deriving the 
Mass-Radius relationship for these objects. The mass (more precisely 
M/R) could be derived from photon redshifts. 
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et al. (2006), is coupling the spin period and field decay as de- 
scribed in equations (|9] - [T2] |. The spin-up torque from a de- 
creasing moment of inertia occurs slowly during the magnetic 
field decay timescale over time ar, and is negligible compared 
to the spin-down due to vortex expulsion. 

The shell eventually reaches an inner radius where it is un- 
stable in the magnetic field geometry, and pieces start to break 
off above Ob- As the pieces fall in along field lines to the polar 
regions and collide with the star, angular momentum is trans- 
ferred. This transfer and the change in moment of inertia, Al, 
as we show below would correspond to a star spinning-up dur- 
ing bursts (i.e. accretion of chunks jfl The change in moment of 
inertia from the increase in mass and radius of the star is given 

by. 



AI 

T 



5 Am 



3 Mqs 



(33) 



where the accreted matter is converted into CFL matter at the 
star's density. Taking into account the angular momentum, AL, 
from the infalling (from R„i) shell material we get. 



AP 
T 



AI 

T 



AL 

T 



5 Amsh 



3 M< 



QS 



(34) 



where a quark star mass and radius of Mqs = \AMq and 
^QS = 10 km is used, and will be assumed for the remainder of 
this section. Making use of equations (fTlfefTsTl for and Am 
respectively one arrives at. 



AP 
~P 



-1.8xlO-V(«^)fi?5 



3^ 
2 Rl 



(35) 



Using equation dZST l. the equation above can be recast into 



AP 
T 



-4.75 X \0-\f{a)PmP-i 



(36) 



where Pio is the spin period in units of 10 s, P-n is the period 
derivative in units of 10"", and we had assumed that R^ » 

If the accretion occurs during a time interval, Afacr, the spin- 
up rate. Pace. - AP/Afacc, can be estimated by dividing both 
sides of equation ( |35] | by Afacc- By noting that Po\i - B^/(3kP) 
we arrive at 



old 



-132/ (a) 



1 hr 



10 • 



(37) 



6.2. Outward Shell Rebound: Spin-Down 

Following the accretion events, during the fast shell rebound, 
angular momentum is transfered from the star to the shell to 
keep it co-rotating, resulting in an increase in the spin-down 
rate of the star. The magnetic coupling between the star and 
the shell implies that the star will lose angular momentum and 
spin-down. As the star moves outward from radius Ri^ back to 



^ In the cases where the shell moves in and out without pieces 
breaking off, there would be associated P variations without SGR 
bursts. 



Rm, conservation of angular momentum of the star/shell system 

implies Lgs.m + Lshell.in = iQS,m + ishell.m- That is 



AP _ 5 Wsh 
P ~^2Mqs 



p2 ^ 

II 



or by using equations ([TI&fTS]) one gets, 
^-+4.1xlO-V(a)PioP-n, 



(38) 



(39) 



which as it turns out is insignificant as compared to mass loss 
from the shell's atmosphere. 

6.3. Shell atmosphere loss: Spin-down 

Observationally, an increase in spin-down has been seen to last 
for ~ 18 days, as in the case of AXP1E2259, and ~ 80 days 
for SGR1900-I-14. In our model this is readily explained by the 
heating of the degenerate shell by the burst, and the subsequent 
release of a portion of it's atmosphere. Because a thin layer on 
the shell's surface is in the non-degenerate phase, a fraction of it 
is blown away via the propeller mechanism, removing angular 
momentum from the system over time. 

To determine the amount of atmosphere present on the shell 
during the quiescent phase, we consider the critical density, 
Psh.nd, of the shell to ensure degeneracy which is found by set- 
ting Psh.eff. = T^Fermi - 76.2 MeVp^|(j)^j j^, wherc the density is 
written in units of 10'° g/cc. So the gas will be non-degenerate 
below densities of. 



Psh.nd - 3.6 X 10 g/cc 



sh.eff. 



1 keV 



3/2 



(40) 



The scale height of the non-degenerate, upper layers, of the 
shell can be estimated to be 



J^I±f^ ^ OA cm RiJ^], (41) 



where g = GMqs/Rf^ ~ 1.87 x 10'7Pf^ cm s"^ is the effec- 
tive gravity. Thus, the mass of the non-degenerate portion of 
the shell is. 



'«sh,nd - 47rPn, sin 6*6 Hsh,ndPsh,n 

~ 3.6xlO'5gm.Pl io( 



,eif. 



IkeV 



5/2 



(42) 



where 0^ - 60° was used. 

Following a bright burst, with luminosity Lb,44 (in units of 
lO'*'* erg s"'), the shell gets reheated to temperatures of the or- 
der Tshit) - 50 keV(Lb,44(f)/^m,io)'^'*- The corresponding at- 
mospheric (non-degenerate) mass in terms of the burst energy 
(in units of 10"*"* erg), is then 



'«sh,nd(f) - 6.32 X 10'^ gm. R^2[l l'^'^^H) 



(43) 



So the trigger mechanism is that the temperature needs to 
be high enough for the atmosphere to leak away from the shell, 
most likely in the form of a pressure driven wind. Once it is 
gravitationally unbound it is kept in co-rotation by the magnetic 
field out to the light cylinder, providing an efficient propeller. 
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The angular momentum per unit mass lost at the light cylinder 
is c^/Q which gives a new spin-down rate for the quark star of, 



-P-ll,new - 2mi2f'io , 



(44) 



where the estimated mass-loss rate, mi2, is in units of 
lO'^gs-i. 

The change in period due to mass lost at the light cylinder 
is also written as. 



AP 



5 Amsh.nd 
2 



M, 



QS 



^QS 



(45) 



where R\c is the radius to the light cylinder Assuming Amsh,nd = 
^'«sh,nd where ( is the portion of the non-degenerate atmosphere 
lost to the light cylinder and Wsh.nd is the peak of the atmo- 
spheric mass during the burst. We get. 



AP 

— = 1.3x lO^'^^L; 
P ^ ' 



5/8 



P,aR 



27/8 



B,44^ 10"m,10 ■ 



(46) 



where in the equation above we use the peak burst luminosity 
which we represent by the subscript "B" as to differentiate it 
from Lb,44(0 ■ 

Because the shell cools, there is decreasing amounts of mat- 
ter available to drive the increased spin-down rate. Thus this 
new rate is temporary and the timescale for it to decay back to 
the previous rate can be estimated by, tp = AOTsh,nd/'«, or. 



Ur -ll.new 



One should keep in mind that in reality the mass loss rate is de- 
creasing in time. Determining better estimates for this requires 
knowledge of ( which is beyond the scope of this paper. 

7. Application to SGRs and AXPs 

7.1. SGR 1806-20: The boxing day SGR 

The most recent and brightest giant flare came from SGR 
1806-20 on Dec. 27, 2004. This flare lasted about 5 minutes 
and had a peak luminosity of about 2 x 10'*^ erg s"'. For a 
distance to SGR 1806-20 of 15 kpc, it is estimated that an 
(isotropic equivalent) energy release of 2 x lO''^ erg occured in 
the spike, and 5 x 10^"^ erg in the tail. 

In our model the energy is readily accounted for The ob- 
served P and P imply (from Eq. |25j that B ~ 1.3 x 10'^ G, 
which translates to an energy of 3.6 X lO'*^ erg for or = 1 and 

= 20 km (see Eq|20]i. 

In the year leading up to the SGR 1806-20 giant flare, 
well-sampled X-ray monitoring observations of the source with 
the Rossi X-ray Timing Explorer (RXTE) indicated that it 
was also entering a very active phase, emitting more frequent 
and intense bursts as well as showing enhanced persistent X- 
ray emission which was, a prelude to the unprecedented giant 
flare. In our model, this prelude would correspond to increas- 
ing amounts of shell pieces breaking off" above 6^ as the shell 
moves in rapidly closer to the star. This situation is likely only 
when the quark star's magnetic field is strong such that it is de- 
caying quickly, implying the object is young, which is the case 
with SGR 1806-20. Furthermore, because ~ 20% of the shefl's 



mass was lost in this giant flare, the leftover smaller shell, once 
re-adjusted to following the event, will have less capability 
to produce flares of the same magnitude as the previous one 
(see eq.(l2ni). 

The sharp initial rise of the main spike in the boxing day 
flare was of the order of 1 millisecond (Palmer et al. 2005). In 
our model, a lower limit on the rise time. 



ARr, 



Vff 



(48) 



is determined by the time it takes a typical chunk of radius 
ARm to cross the star surface, and be converted to CFL quark 
matter. Larger rise times would correspond to groups of shell 
pieces accreting simultaneously. To a first approximation we 
assume that the shell piece crosses the quark star surface at 
maximum (free-fall) velocity of, Vfj- = -y/GMgs /Rqs- Making 
use of equations (I2ll20l&l25]l we get. 



frise.g > 0.5 flS 



VPiof_ii(?7o.i/(g)) 

f1/4 
^b,47 



1/4 



(49) 



where the burst energy, given in units of 10 erg, is represen- 
tative of a giant burst or the case where a > 1 (denoted by 
subscript "g" in the equation). 

For regular bursts (a « 1 ) we find 



(47) frise,r > 0.9 flS 



Vf IpP-ll (/;o.lQ'-6) 
^b,40 



1/4 



(50) 



where a^d — a/ 10 was used and the burst energy is given in 
units of 10"^" erg. However, in the case of small bursts the total 
rise time can be larger, as multiple shell pieces are likely to be 
falling into the star simultaneously. 

7.2. SGR 1900+14 

The rotation period has been studied in detail for SGR 1900h- 14 
during bursting (i.e. Woods et al. 1999 and Palmer 2002), 
which allows us to test these features of our model. The in- 
crease in period as given in our model (Eq. |46] | following the 
burst is, AP/P ~ 3.9 x 10"^ x ^^^''(q using observed values 
for energy and period. Also, observed values for period and the 
new period derivative following the August 27th burst give a 
mass loss rate of m ~ 5 x lO'-' g/cc, and an upper limit for the 
recovery time of tp ~ 15 days x ^^^(y- 

The magnetic equilibrium radius was estimated to be 80 km 
due to the presence of iron line emission (Strohmayer & 
Imbrahim 2000) in which case for ( ~ 0.004 we get AP/P ~ 
1.5x10 "* andfp ~ 67 days. However as discussed in § 17.2. Fl for 
an ionized atmosphere, redshift correction could move the ra- 
dius to 30 km in which case a ^ ~ 0.1 gives AP/P ~ 1.4 x 10"'* 
and fp ~ 61 days. This is in good agreement with what can be 
inferred from observations. 



7.2.1. Iron line emission during bursts 

Observations of Fe fluorescence lines in SGR 1900H-14 
(Strohmayer & Ibrahim 2000) are a key feature in our quark 
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star model. As discussed in Ouyed et al. (2002) and Keranen, 
Ouyed, & Jaikumar (2005), Fe and heavier metal production 
in the ejecta is significant during the formation of a quark star. 
Strohmayer & Ibrahim state that the distance from the star's 
surface, where the emission is created needs to be at least, 
h ~ 80 km, to account for the lack of redshifj^. This dis- 
tance gives an excellent confirmation to our model, as it cor- 
responds to the distance of the shell which possess an Fe-rich 
atmosphere. 

The atmosphere's column density can be calculated from 
equation (l42ll to be. 



A^Fe = 3.8 X lO^"* cmr^ 



m,10 



eff. 



IkeV 



5/2 



(51) 



This is sufficient to produce a strong Ka iron emission line 
(equivalent width > 100 eV), when illuminated. Thus in our 
model, the Fe flourescent line is produced during bursts by the 
illumination of the shell's atmosphere. 

7.3. X-ray to Infrared Flux Ratio 

Another observed feature in both SGRs and AXPs is the X-ray 
to infrared flux ratio correlations. The X-ray flux during the 
quiescent phase (induced by the magnetic field decay) and/or 
during the bursting phase (from accretion events) can be ther- 
mally reprocessed by the shell into the far infrared. However, 
the efficiency is too low due to the small shell size and cannot 
account for the observed values of (vFy)x/(vi^v)iR which are 
~ 150 for SGRs and ~ 1500 for AXPs (e.g. Israel et al. 2003). 

Alternatively one expects a much higher infrared flux re- 
lated to vortex annihilation, which is always present in this 
model. The mechanism is synchrotron emission from mildly 
relativistic electrons accelerated by the magnetic reconnec- 
tion events induced by vortex annihilation (see simulations in 
Ouyed et al. 2005). The X-rays are from the high energy elec- 
trons and the infrared from the low-energy electrons. 

For a power-law electron energy distribution N{E) oc 
the optically thin synchrotron flux is given by Fy oc y-^r-^y^ 
(Longair 1992). This implies 



(vFy)lR VIR 

For a typical index y = 2 this yields 



{yFy)u 



140. 



(52) 



(53) 



In this model both Lx and Ljr are proportional to (see 
Eg. |30]l which also implies Li^lLx is a constant. Synchrotron 
self-absorption at lower frequencies, low-energy electron 
losses, or a smaller index y, all increase the Fx/ Fir ratio, and 
may be the reason why AXPs show larger ratios than SGRs. 
During bursting episodes the vortex annihilation rate increases 
due to the spin-up induced by accretion increasing both Lx and 
LiR, thus they remain correlated during bursts. 



* If the line is from ionized iron, the rest energy can increase up to 
6.7 keV thus a redshift up to 5% is allowed, reducing the lower limit 
on the distance of the iron emitting gas from the star to ~ 30 km. 



7.4. The hard spectrum/component in our model 

Three of the four known SGRs have had hard spectrum (with 
photons in the MeV energy) giant flares. Before showing how 
this can be accounted for in our model, recall that in the case 
of quark stars the surface emissivity of photons with ener- 
gies below hcjp ^23 MeV (Wpi electromagnetic plasma fre- 
quency) is strongly suppressed (Alcock, Farhi, & Olinto 1986; 
Chmaj, Haensel, & Slominsli 1991; Usov 1997). In Vogt, Rapp, 
& Ouyed (2004) it was shown that average photon energies 
in quark stars in the CFL phase at temperatures T are ~3T. 
Therefore, as soon as the surface temperature of the star cools 
below Ta - hajp/3 ^7.7 MeV, the photon emissivity is highly 
attenuated. This is studied in more details in Ouyed, Rapp, 
& Vogt (2005) where it was demonstrated that for tempera- 
tures above 7.7 MeV, neutrino cooling is dwarfed by the pho- 
tons; i.e., photon emission/cooling dominates as long as the star 
cools from its initial temperature Tq > 7.7 MeV to = 7.7 
MeV. For temperatures below 7.7 MeV, cooling is dictated by 
the slower neutrino processes. 

To a first approximation, the increase in the star's tempera- 
ture, AT^Qs, following the accretion of shell pieces can be writ- 
ten as. 



47iRq^AR(1 - sin0B)cvArQs - r]AmshC^ , 



(54) 



where Cv = 7.8 x 10 I'^j.y erg cm" K" is the star's spe- 
cific heat as given in Ouyed, Rapp, & Vogt (2005). Also, AR 
is the thicknes^ of the heated region in the polar caps. Since 
the pieces fall off above the neutral line (i.e. > Ob), the polar 
regions are the ones that are primarily heated in the process. 
The corresponding increase in temperature (assuming the ini- 
tial star temperature is in the keV range) is then. 



Tqs ^ 10 MeV 



j.3/16 
"b,45 



((l-sin0B)C/'^'(«)4s,io^i5) 



1/4 



(55) 



For Ob - 60°, a - I, and B15 ~ 1, this implies that only ac- 
cretion events with energy exceeding Eb.c ~ 2.8 x 10"^^ erg are 
capable of reheating the star above 7.7 MeV, allowing MeV 
photons to be generated and escape the star. Thus our model 
naturally explains why only the rare giant bursts emit a hard- 
spectrum. 

The escaping photons are thermalized and cool the star at a 
rate given by 47r(l - sm0B)RQgARcydTQ5/dt 
which leads to a cooling time of. 



4nRl,(rT'^, 



fc.QS ~ 0.03 ms 



(l-sin0B)/7,y,V^V)fii5. r, 



I OA 

7I/4 
"b,45 



, -QS 
log — ■ 

^ a 



(56) 



As the high energy photons escape into the magnetosphere they 
interact with lower energy photons and (e^e ) pairs, and so be- 
come thermalized to about ~ 1 MeV. The (e^e") pairs could be 
due to the radiation from the magnetic field decay which oc- 
curs outside the surface of the quark star. We note that the radi- 
ation produced by accretion is from the conversion of baryons 



' The heat penetration depth is given as AR = cjg where Cs ~ c/ V3 
is the CFL sound speed and tff is the free fall timescale given by equa- 
tion We get AR ^ 10* cm Bis x (?7o.i/(a)/£b,47)"^ 



Ouyed et al.: Quark-Nova remnants: application to SGRs, AXPs and XDlNs 



11 



to quarks at the surface which also releases photons at the sur- 
face (half of this radiation heats the quark matter inside the star 
while the other half is radiated promptly). We have not stud- 
ied these emission mechanisms in details but we are currently 
investigating mechanisms for thermalizing the escaping radia- 
tion, and the nature of the resulting spectrum. We can simply 
argue for now that the radiation mechanisms should be simi- 
lar to those for a pair-dominated fireball in an optically thick 
environment. 

7.5. Summary of outstanding SGR/AXP questions 

We summarize this section by attempting to test our model 
against the open issues related to SGRs/AXPs as discussed in 
the literature (i.e. Israel 2006). These open issues are enumer- 
ated below. 

1 . Are the SGRs/AXPs engines born with milliseconds peri- 
ods? 

Allen & Hovarth (2004) and Vink & Kuiper (2006) give 
good evidence in two cases for normal energy super- 
nova shells ai-ound an AXP (1E1841-045) and an SGR 
(0526-66), which implies birth periods larger than tens of 
milliseconds. In our model there is no need for rapid rota- 
tion at birth. 

2. What diff'erentiates the two types (giant and regular) of 
bursts? How does the boxing day, SGR 1806-20 event fit 
in this picture? 

A giant burst is due to the shell losing a larger fraction 
of its mass as it moves towards the star, parameterized in 
our model by a ~ 1. Regular bursts, a <s 1, are due to 
smaller pieces breaking off the neutral line as the shell os- 
cillates around its equilibrium position. For SGR 1806-20 
we argue the boxing day event is one of the first events ex- 
perienced by this object following its birth. Moreover, our 
model predicts (via eqn.l2TI) that we should continue to see 
a decreasing trend in both the burst intensity and frequency 
going from younger to older objects. 

3. Why are the IR and X-ray variability correlated during 
flares? Is the far-IR emission due to a passive disk? Why 
passive? Why disk? 

The far-IR emission and the correlation with the X-ray 
emission can simply be explained as synchrotron emission 
from the high-energy and low-energy electrons (see § 17. 3b . 
Our disk is in fact an iron-rich shell. It is passive because 
it is degenerate and for most of its lifetime it remains in 
equilibrium at R^^, co-rotating with the star 

4. Is there any connection between AXPs/SGRs and high-B 
radio pulsars? 

AXPs/SGRs are quark stars and high-B radio pulsars are 
neutron stars that have not gone through a QN phase, prob- 
ably because their core densities never reached deconfine- 
ment values (see Staff, Ouyed, & Jaikumar 2006). 

5. What is the origin of the hard X-ray spectrum in 
AXPs/SGRs? 

The hard X-ray spectrum in our model can be explained 
as MeV photons generated in the outer layers of the star 
These photon bursts can only occur for accretion events ca- 



pable of heating the star above 7.7 MeV. We find that only 
bursts with energies above lO"*"* erg can do so (see § 17.4b . 
thus explaining why only the giant SGR flares show a hard 
spectrum. 

6. Why do the spin-down ages and the supernova ages differ? 
In our model the difference can be explained by the time 
it takes the neutron star to reach deconfinement/nucleation 
densities as discussed in § 14.41 Simply put, the supernova 
age is the time for the neutron star to reach quark deconfin- 
ment densities and experience a quark-nova, plus the time 
needed for strange quark nucleation (fsN = ^qn + fnuci.)- 

7. What are the progenitors of AXPs/SGRs and why are they 
located in dense ISM? 

Massive stars (near the black hole line) exploding in high 
density ISM are most likely to lead to massive compact 
remnants. The high density ISM will confine the massive 
progenitor winds much closer to the star, causing the decel- 
eration of the blast wave, and initiating the reverse shock 
inside the remnant (Truelove & McKee 1999). This would 
lead to more massive compact stars which are more likely 
to turn directly into quark stars. 

8. XDINs in our model 

In our model X-ray Dim Isolated Neutron Stars (XDINs) are 
old SGRs/AXPs that have gone through their most active burst- 
ing phase and are left with a thin shell in stable equilibrium at 
/?ni- We start by summarizing the observed and measured fea- 
tures of XDINs before we apply our model to these intriguing 
objects. 

8.1. Properties of XDINs 

These dim (Lx ~ 10"" erg s"') isolated neutron stars are nearby 
at around 100- 300 pc and show no SNR association. Three 
of them have known proper motions that are too fast to ac- 
crete. The main common properties of the "magnificent seven" 
can be summarized as follows (see Haberl 2005 and references 
therein): (i) The blackbody: The X-ray spectra of XDINs ob- 
tained by the ROSAT PSPC are all consistent with blackbody 
emission. The soft X-ray spectra have a temperature in the 
range of 40-110 eV. They show no non-thermal component; 
(ii) The Absorption lines: The XMM-Newton spectra can be 
best modeled with a Planckian continuum including a broad, 
Gaussian shaped absorption line. The line centroid energies are 
in the range 100-700 eV. The depth of the absorption line (or 
the equivalent width) was found to vary with pulse phase. It 
has been suggested in the literature that these absorption lines 
can be best explained as proton cyclotron resonance absorption 
features in the 0.1-1 keV band with field strength in the range 
of 2 X 1013-2 X lO'"* G (Zane et al.2001; Zavlin&Pavlov 2002) 
with the line broadening explained as due to the variation of 
the magnetic field over the neutron star surface; (iii) The op- 
tical excess: They show optical excess compared to the X-ray 
blackbody. In other words at optical wavelengths they show 
a factor of about 3-14 excess when compared to the extrap- 
olation from X-rays (Pons et al. 2002; Motch&Haberl 1998; 
Haberl et al. 2004; van Kerkwijk et al. 2004); (iv) The Lack 
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of radio-emission: It has been argued that the lack of pulsed 
radio-emission is because their radio beam is very naiTow due 
to the large light cylinder radius (i.e., large periods). However, 
there exist radio-pulsars with similar magnetic field strengths 
and periods (e.g. Camilo et al. 2000) that are active in the ra- 
dio. In our model, the lack of radio pulsation is due to the fact 
that quark stars in the CFL phase, unlike neutron stars, become 
aligned-rotators following the QN (see Ouyed et al. 2005). We 
now go on to discuss the remaining properties. 

8.2. The two-component blackbody 

As discussed in ^5.11 the first blackbody temperature is set by 
the X-ray continuum luminosity from the surface of the quark 
star, 



Tbb ^41.0 eV 



QS.IO 



1/2 



(57) 



Note that the spin-down is in units of lO"'-' s/s reflective 
of what has been measured for XDINs. The observer would 
see an emitting region with a corresponding blackbody radius 
^BB = ^QS VI - sin^B ~ 4.4 km using our fiducial value 
0s = 60°. Interestingly, these are the same values as inferred for 
RX J1856 (e.g. Burwitz et al. 2003) and RXJ0720 (e.g. Haberl 
2004). 

Similarly, the effective shell temperature (the second black- 
body) is 



sh.eif . 



41.0 eV 



,10 



1/2 



QS.IO 



Rn 



1/2 



(58) 



8.3. Absorption lines 

Since R^ > Rqs in our model it implies Tsh.eiF. < T'bb which 
means that the shell will act as an absorber of the hotter X- 
ray blackbody. The maximum density, ps.nd, of the shell below 
which the gas is non-degenerate is found by setting Tsh.eff. = 
T^Fei-mi, or, psh,nd - 54 g cm"^ X (rsh,eff./100 cVf^^. The scale 
height of the non-degenerate, upper layers, of the shell can be 
estimated to be 



i/.,„a~^^~0.01cm7^^,(^^ 



(59) 



I^Ffiinug, \100eV/ ' 

implying a mass of the non-degenerate portion of the shell of 

mshM = ^^Rm sin H^hMPsM (60) 
with an optically thick column density. 



N^, ^ ^r'^^-"" -2.1x10"^ cm-^/;^ 



T'sh.eff 

(56mH) ' "'"'llOOeV 

Tbb 



5/2 



(61) 



2.1 X 



inl8 -2 d3/4 / BB \ 

10 cm R ' ,„ 

"•'"UOOeV/ 



Iron photospheric models have been calculated by Rajagopal et 
al. (1997) for magnetic field strengths as derived for XDINSs 



in our model (i.e. B - yfJiiPP ~ 10'^ G). Despite the fact 
that these calculations are only done in Hartree-Fock approx- 
imation the results show absorption features in the blackbody 
spectra, at energies 300 eV and above, reminiscent of the lines 
in XDINSs (see Figure 3 in Rajagopal et al. 1997; see also dis- 
cussion in Neuhauser et al. 1987). Careful attention to Figure 
3 in Rajagopal et al. (1997) shows that the lines start to dis- 
appear at temperature between 30 eV and 80 eV. Interestingly, 
this temperature efi'ect may explain the absence of absorption 
line in the two XDINSs showing the lowest blackbody tem- 
perature namely, RX J0420.0-5022 with kT^,^ ~ 40 eV and 
RX J1856.5-3754 (the brightest or closest of the XDINs) with 
^T^bb ~ 60 eV. We thus suggest that the iron shell in our model 
provids the conditions necessary to explain the absorption lines 
as iron lines. In this respect, the iron atmosphere explanation of 
the spectrum suggests that younger (i.e. hotter) objects of the 
same type should be observed with a richer line spectrum. 

8.4. The optical excess 

The shell's solid angle, 4?: sin 0b, is essentially providing the 
excess optical emission. In the simplest case, one assumes a 
scattering atmosphere isotropizes the optical emitted flux from 
the QS (i.e. blackbody tail) and the shell. Then the ratio be- 
tween the total optical flux and the optical contribution from 
the tail of the blackbody can be expressed as 



opt., tot. 
^ opt., tail 



^ 1 + 



opt.,sh 



= 1 -t- 



-^opt.,tail 

An sin 0^ 



(62) 



1 



47r(l - sin ^b) 1 - sin 0b 



In other words, the optical excess is a direct measure of the an- 
gle subtended by the shell. This also means a one-to-one corre- 
spondence between the optical excess and the size of the emit- 
ting spot in our model expressed as 



opt., tot. 



^ 1 _ (Rqs 

^ opt., tail l-sin^B \Rbe 



(63) 



For example for RXJ1856 R^b ~ 4.4 km has been derived 
from observations. Assuming a quark star radiu^ of about 
10 km, our model would then predict an optical excess of 
(10/4.4)^ ~ 5 which is very close to the factor 5-7 measured 
from the current optical data and the LETG spectrum (Haberl 
2004). The radius of the emitting area in the case of RXJ0720.4 
has also been estimated to be 4.4-4.8 km (for a distance of 300 
pc) which should give an optical excess of ~ 5 in our model 
and similar to what is observed (see Table [TSIi. 



* It was suggested that the minimum radius of RXJ1856 might ex- 
ceed 14 1cm thus favoring stiff equations of state (Triimper 2005). We 
argue that the inferred radius is in fact the location of the iron shell. 
Indeed, the temperature of the cool component was measured to be 
< 33 eV at the 3cr level (Burwitz et al. 2003). In our model, it implies 
r,h.eff. < 33 eV 01 R,, > (rBB/r,h.eff.)'i?Qs - (60./33.)2ijQs. That is, 
> 40 km assuming Rqs ~ 10 km. 
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9. Conclusion 

In this paper we present a new model for SGRs and AXPs with 
possible appUcations to XDlNs. This novel idea relies on the 
formation of bare CFL quark stars (within the quark-nova sce- 
nario) as the underlying engine, with the parent neutron star's 
crust material surrounding it. Despite the simplifications we 
have employed to study the dynamics of this ejected shell, 
specifically how it reaches the equilibrium radius and how it 
evolves once it is there, our model has a number of attractive 
features that can account for many observed SGR, AXP, and 
perhaps XDIN properties. 

Missing from our study is the proper treatment of three- 
dimensional instabilities acting on the shell, such as Raleigh- 
Taylor or oscillations. Due to the high conductivity of the 
shell we expect the Raliegh-Taylor to not operate on very fast 
timescales, and any radial oscillations to be severly damped as 
the conducting shell will maintain the magnetic flux enclosed 
within. However, azimuthal oscillations may have a significant 
effect. For a more accuarate description, the dynamics of the 
shell will require numerical simulations, and is left as future 
work. 
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Table .1: Observed features of SGRs and AXPs (Only AXPs and SGRs 
with measured P and P are listed)^ 



name 


agesNR (kyr) 


P(s) 


P-u (s/s) 


-I- 


ix(ergs-') 


^g.burst (erg) 


fg,rise (ms) 


^n.burst (erg) 


fn,rise (ms) 


AP/P 


grl 806-20 


17 ± 13 


7.49 


8.3 


1.28el5 


3.6 X lO^*- 


~ 3.7 X 10'*'' 


< 1 


10^0-104' 


- 


- 


,grl900+14 


20 ± 10 


5.17 


6.1 


9.12el4 


2 X 10^'' 


> 1.2x 1044 


<4 


lO-^C-104' 


- 


+ 10-4" 


,gr0525-66 


10 ±5 


8.05 


6.5 


1.18el5 


8 X 10'' 


5.2 X 1044 


< 2 


10-*"- 104' 


- 


- 


1E1048-5937 


20 ± 10 


6.45 


3.81 


8.05el4 


3.4 X lO-** 


- 


- 


2.7 X 10472.8 X 104' 


21/5.9 


- 


iXteJ1809-1943 




5.54 


2.06 


5.49el4 


1.6 X lO^*- 












)lE1841-045 


1.5 ± 1 


11.77 


4.13 


1.13el5 


2.3 X lO-'^' 












:pJO 100-72 11 




8.02 


1.88 


6.31el4 


1.5 X lO-*^' 












IXS 1709-4009 


17 ± 13 


11.00 


1.94 


7.51el4 


6.8 X 10^^ 












)4U0142+615 


no SNR 


8.69 


0.196 


2.12el4 


3.3 X 10^4 












)lE2259+586 


10 ±7 


6.98 


0.0484 


9.44el3 


3.8 X lO-*"* 






[3.4 X 10-"]" 


(1-100) 


-4 X 10-*'"'" 



im the SGR/AXP online catalogue at : http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 
ejects are listed in order of increasing estimated age (using our model; see § l4.4l l: 
isient AXPs; 

;rease in spin-down rate with /"new ~ 2.3/'oid for about 80 days; 

y sudden spin-up during the burst with /"new ~ 2/'oid for about 18 days; 

ived in our model from B - \3kPP; 

bursts with energies ranging from 3 x lO-'" to 5 x lO-'^. 



Table .2: Observed features of the X-ray dim radio-quiet isolated neutron 
stars** 



ID 


period (s) 


^-13 (S/S) 




Lv (erg s"') 


/vTbb (eV) 




£cen. (eV) 


Oplical Excess 


y^bb (km) 


RX J0420.0-5022 


3.45 


< 92 


- 


2.7 X lO-*" 


44 


- 


- 


< 12 


- 


RXJ0720.4-3125 


8.39 


0.69 


3.93el3 


2.6 X 10^' 


85-95 


40 


280 


5 


4.4-4.8 


RX J0806.4-4123 


11.37 


< 18 




5.7 X 10^" 


96 




306-430 






RBS 1223 


10.31 


1.12 


5.52el3 


5.1 X 10* 


86 


150 


230-300 


<5 




RXJ1605.3+3249 










96 


120 


450-480 


< 14 




RXJ1856.5-3754 










60 






5-7 


4.4 


RBS 1774 


9.44 








101 




~ 700 







See Haberl (2005) and references therein. 



